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Abstract 



The question of the conditions under which ID systems support 
extended electronic eigenstates is addressed in a very general con- 
text. Using real space renormalisation group arguments we discuss 
the precise criteria for determining the entire spertrum of extended 
eigenstates and the corresponding eigenfunctions in disordered as well 
as quasiperiodic systems. For purposes of illustration we calculate a 
few selected eigenvalues and the corresponding extended eigenfunc- 
tions for the quasiperiodic copper-mean chain. So far, for the infinite 
copper-mean chain, only a single energy has been numerically shown 
to support an extended eigenstate [ You et al. (1991)] : we show ana- 
lytically that there is in fact an infinite number of extended eigenstates 
in this lattice which form fragmented minibands. 



In recent years extended electron states have been found in a variety 
of diverse ID systems [1-14] ranging from disordered systems on the one 
hand to various kinds of quasiperiodic systems on the other. In most cases, 
such states are identified through extensive numerical computations [1-8] . A 
general analytical understanding of this problem is still lacking [9-14]. In 
this paper we like to address this problem and we wish to see under what 
physical conditions extended states are expected to occur in disordered or 
quasiperiodic systems. 

Consider an elementary example - an Anderson localisation problem in- 
volving the random distribution of A atoms on a host B-atom lattice. It is 
well-known |3], |l(J that there is an extended eigenstate at a certain energy in 
this system if the A-atoms always occur in pairs, i.e., in the form of dimers. 
This is an instance of a definite kind of correlation between atoms leading 
to extended states in a disordered system. In this case, the correlation con- 
sists in the fact that an A atom always occurs as a member of a dimer. We 
may generalise this idea as follows. Suppose that there is a certain cluster 
of atoms which is distributed in some manner on a linear chain. Let this 
cluster consist of a finite repitition of a certain subcluster of atoms. Thus in 
the above example, a cluster (an AA pair) is made up by repeating a sub- 
cluster (an A atom) a finite number (two) of times. If this subcluster is to be 
found only within the clusters and nowhere else in the chain, then it can be 
shown [0 that for certain discrete energy values these clusters offer identity 



contributions to the total transfer matrix [fig] . At these special energies we 
may disregard the presence of the clusters themselves, and if it happens that 
the remaining part of the lattice is a periodic one, then evidently this lattice 
will support extended states at these energies. 

The special energies which support extended states are a property of 
the cluster only, and do not depend on the lattice as a whole. Thus, in 
our example of randomly distributed A-type dimers on a B-type host lattice 
(RDL), the special energy value E = e& at which there is an extended state 



is determined from the consideration of the AA cluster alone [0. At this 
energy value these AA dimers make identity contributions to the total trans- 
fer matrix and the chain effectively consists only of B atoms. We now wish 
to emphasize that many other ID lattices containing these AA dimers may 
have extended states at this same energy. We may illustrate this by consider- 
ing an ordered chain consisting only of A atoms and comparing its spectrum 
with that of the RDL. In the ordered chain the AA correlation is trivially 
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present, and therefore this system supports an extended state at the same 
energy. However, and this the crucial point we wish to make, while all eigen- 
states at other energies are localised ones in the RDL, in the ordered chain 
of A atoms we have a whole band of extended states, the existence of which 
cannot be inferred from a consideration of the AA dimer alone. Although 
this band can be found directly by using Bloch's theorem, we do encounter 
many systems also possesing bands or minibands of extended states [5-14], 
where this theorem cannot be used as there is no periodicity. An example of 
this is the quasiperiodic copper-mean chain (CMC). 

The clusters contain only part of the information about the entire lattice, 
and therefore, provide a common set of energy levels for a variety of systems 
which contain these clusters but which otherwise differ among themselves in 
their long range behaviour. In general, these systems possess correlations 
beyond the level of what is contained in the clusters themselves, which are 
responsible for the detailed differences in electron spectra among them. In 
order to bring out these differences we have to examine correlations among 
larger and larger blocks of atoms in the chain, by looking at the system at 
larger and larger length scales. In this process the long range features of the 
lattice, which may include periodicity, or may be self-similarity, get gradually 
included in the determination of the electronic spectra, and we end up by 
obtaining the entire spectrum of extended states. 

In this paper we present an approach which incorporates this point of 
view for finding the full spectrum of extended states in ID lattices, periodic, 
quasiperiodic or disordered. The most natural method of including the effects 
of correlations at all length scales is the real space renormalisation group 
(RSRG) method, which is adopted here. 

Although the discussion so far has been perfectly general, from now on we 
illustrate our ideas by focussing attention on the CMC. A portion of a CMC 
is shown in Fig. 1. The sequence in which long (L) and short (S) bonds are 
arranged in a CMC may be obtained from the substitution rule L — > LSS 
and S — > L starting with, say, L. For decribing the electron states in this 
lattice we use the tight binding hamiltonian 

# = E e *M*l+E^l i >01 C 1 ) 

i <ij> 

As shown in Ref. J7|, for implementing the RSRG method in CMC we 
have to identify four types of sites, namely, a, /3, 7 and 5 corresponding 
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respectively to the L — L, L — S, S — L and S — S vertices (Fig. 1). The site 
energies in Eqn.(|l]) assume the values s a , Ep, e 1 and e§, and there are two 
different values of the nearest neighbour hopping integrals, t L and t S - 

The eigenf unctions for such an ID lattice may be calculated by the stan- 



dard transfer-matrix method [15fl . The discretized Schrodinger equation can 
be cast in the form 4> n+ i = M n (p n , where 




E-en 



and M„ = tn,n+i t n , n+1 

\ 1 J 

Here ip n denotes the amplitude of the wavefunction at the n-th site and M n 
is a 2 x 2 transfer matrix. In the CMC, <p n is related to <fio by a product 
of four types of transfer matrices M a , Mp, M 7 and Ms following the copper 
mean sequence. 

By inspecting the CMC (Fig. 1) we see that the a-sites always occur in 
pairs, whereas, ft, 7 and 5 always form a triplet ftSj. Thus, for this lattice, 
the pair of sites aa constitutes a cluster. If we now renormalise the lattice by 
applying the deflation rule @ LSS — > L and L — > S, we get a scaled version 
of the original lattice with scale factor a = 2. In the renormalised chain we 
again find a sites occuring as aa pairs, and ft, 7 and S sites forming ftSj 
triplets. This means that the aa clustering effects are also present on this 
inflated length scale, and therefore, by induction, present at all length scales. 
It is important to appreciate that the aa clustering at larger length scales 
amounts to subsuming, so to say, the effects of larger and larger segments of 
the original chain into the a subclusters, thereby allowing larger and larger 
portions of the chain to contribute to the special energies for the extended 
states. 

It is clear from the above discussion that the aa correlation at all length 
scales can only be revealed by renormalisation group methods. While in the 
original lattice, this correlation is directly visible (Fig. 1), higher order cor- 
relations due to renormalised aa pairs imply underlying complex correlation 
between atoms, which is not apparent from mere inspection of the original 
lattice. 

In a copper-mean lattice the string of transfer matrices typically looks 
like, 

■ • • Mp 5y M 2 a M^M 2 a Mp S , ■ • • , (2) 
where, Mp &1 = M^M h Mp- 
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We notice that the matrix M a is unimodular, and so we can apply a well- 
known result, due originally to Cayley and Hamilton, for the m-th power of 
a 2 x 2 unimodular matrix M : 

M m = U m ^(x)M-U m _ 2 (x)I (3) 

where, x = (l/2)TrM. U m (x) = sm(m+l)9/ sin 9, with 9 = cos -1 x is the m- 
th Chebyshev polynomial of the second kind. Let the matrix M m correspond 
to the transfer matrix of a cluster composed of m identical subclusters, each 
described by a 2 x 2 unimodular transfer matrix M. For a value of the energy 
E for which U m _i(x) becomes equal to zero, we have M m = — f/ m _ 2 (x)J, 
that is to say, the transfer matrix for this cluster behaves essentially as the 
identity matrix at this energy. Thus at this energy, the lattice does not 
'feel' the presence of the clusters defined by the transfer matrix M m . If the 
remainder of the lattice forms a periodic chain, there will be an extended 
state at this energy provided this energy is an allowed one. For allowed 
states wavef unctions do not diverge at infinity |T5 



For the CMC, therefore, we have to put m = 2 in Eqn.(^|). We see 
that Ui(x) = for E = e a , and therefore, for this energy M 2 = —I as 
Uq(x) = 1. Thus, for this energy, the string of transfer matrices given by 
Eqn. (2) is composed only of matrices Mps^. This effectively describes an 
'ordered' linear chain composed of unit cells each containing three types of 
atoms (3, 5 and 7. Now, if the condition (l/2)|TrM / 35 7 | < 1 is satisfied for 
E = e a , then this energy is an allowed one (see Ref. ||15||), and we shall have 
an extended eigenstate for the whole system. 

To determine all the other extended states, let us now consider all the 
successive renormalised versions of the original chain. The self-similarity 
of the lattice implies that one can apply the considerations discussed so 
far in this paper to every such renormalised version. At any stage i of 
renormalisation, extended states will be found at energies for which E = 
where e^'(E) is the renormalised site-energy of a site at this stage. The 
expression for can be obtained by iterating the recursion relations M 

e' a =e 1 + Q(t 2 L Pp + t 2 s P s ), e' p = e 1 + Qt 2 s P 6 , ^ = e a + Qt 2 L P p , 

e' 5 = e Q , t' L = Qt 2 t L PpP s , t' s = t L} (4) 

where Q = (1 - t^PpPs)" 1 and P» = (E - with i = a, (3, 7, 5. 
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The roots of the polynomial E — for every £ yield a new set of energy 
values for which eigenstates are extended. The number of allowed energy 
values increases with the progress of iteration. The totality of all these ener- 
gies constitutes the entire spectrum of extended states for the CMC. These 
energies form minibands in energy space, the existence of which was only 
detected previously through numerical work . 

The CMC may be divided into three sublattices Oj, f2 2 and f2 3 (Fig. 1), 
each being a scaled version of the original lattice 0. The recursion relations 
in Eqn.(§|) actually represent renormalisation of the original chain to the fli 
sublattice. We can define similar RSRG transformations for the 2 arid Q% 
sublattices. Let us denote the transformations for these three sublattices by 
Ti, T2 and T3. Previously we have determined e^) by applying Ti transfor- 
mation I times successively. In general e^} is to be determined by applying 
a string of I operators consisting of Ti, T 2 and T 3 in any arbitrary sequence. 



A string of T's represents a genealogical path in the sense of Ref.[[U; 
and all the extended eigenstates corresponding to a given string of T's can 
therefore be considered as members of the same 'family'. These states are 
characterised by the same length scale, and they belong to the same point on 
the genealogical tree. Now it may happen that the states belonging to the 
same family appear in different regions of the energy spectrum, that is, they 
may belong to different minibands. In other words, two neighbouring states 
belonging to a small energy interval AE may belong to entirely different 
families. This indicates that the minibands are composed of extended states 
which are, in general, characteristically different from each other, forming 
dense sets of discrete energies. 

We now present some results of our calculations for the CMC. Taking 
the hamiltonian parameters as e a = Ep = e 7 = e$ = 0, and ts/ti = 2, the 
first few stages of renormalisation yield the extended state energy eigenval- 
ues E<® = 0, E^ = ±3, E® = ±1.3637980258 and ± 3.0232523786, and so 
on. Only the case = has been detected earlier M through numerical 
calculations. In Fig. 2 we have displayed wavefunctions for a 20th generation 
copper-mean lattice consisting of 699051 bonds. The first column of Fig. 2 
shows the amplitudes on the first 450 sites of this chain, for three selected 
energy values belonging to three different genealogies. The next column in 
this figure shows the amplitudes on the last 450 sites on this chain, for the 
same three energy values. Examining across a row we notice that the 'local 
pattern' of the wavefunction is the same throughout the chain, with the am- 
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plitudes oscillating in an aperiodic manner within definite bounds without 
decaying, unlike the behaviour displayed by a critical or localised eigenfunc- 
tion. The pattern itself is different for the three different rows. Although 
these states are extended, they are not Bloch-like. For E = Gumbs and 
Ali |L7j found that the amplitudes decay with a power law. However our 



analysis and explicit evaluation of the wavefunction clearly shows that their 
conclusion is incorrect. 



As has been shown in an earlier work ||11|| , in a quasiperiodic lattice the 
flow pattern for the hamiltonian parameters under renormalisation at energy 
eigenvalues for extended states displays a regularity which is absent in the 
non-extended case. In the CMC, this regularity associated with an extended 
state shows up in the site energies for the (3 and 7 sites becoming equal, after 
a certain step of iteration, the equality being maintained in all subsequent 
steps. The level at which the f3 and 7 site energies become equal for the first 
time depends on the energy value. Furthermore, this level is the same for all 
energy eigenvalues belonging to the same family. 

The ideas discussed in this paper are very general. Thus an ordered chain 
of A atoms (e, = e and = t) 

may be regarded as consisting entirely of AA dimers, which remains true 
at all subsequent length scales upon renormalisation of the chain. At any 
level, say n, of iteration the energy eigenvalues are given by the following 
analytical expression 



E {n) =e±ty2± v / 2± v / 2±---n times. (5) 

The eigenvalues E^°\ E^\ E^ 2 \ . . . E^ taken together constitute the entire 
band ranging from e — 2t to e + 2t. The same approach applies to the 
binary ordered chain of A and B atoms. Here again, the band edges and the 
band gap extending from ea to eb due to the level repulsion are correctly 
located. The problem of a random distribution of dimers on a lattice has 
been already discussed, there being a single energy at which there is an 
extended state. This is related to the fact that there is no deeper level of 
correlation in this chain than what is contained through the existence of 
dimers. If a disordered chain possess higher order correlations among the 
sites in the sense discussed in this paper then there exist extended state at 
several energies. This may explain the various observations reported in the 
literature on disordered systems [1-4]. 
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Figure captions : 

Fig.l The aa and fiS'j blocks in the copper-mean chain and also in its sub- 
lattices. 

Fig. 2 The amplitudes of the eigenfunction on the first 450 sites and the last 
450 sites of a 20th generation CMC having 699051 bonds. Here e a = 
Ep = e 1 = e$ = and ts/ti = 2. The 1st, 2nd and 3rd rows correspond 
to E = 0, -1.3637980258 and -2.6255768195 respectively. All energies 
are measured in units of t L . We choose the boundary condition as 
■^0 = and ipi = 1. 
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